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A search for neutral Higgs bosons produced in association with b-quark(s) and 
decaying into a pair of b-quarks is performed with the CMS detector at LHC. 
The Higgs boson signal is expected to emerge as an excess in the mass spectrum 
of two b-tagged jets. Dedicated triggers with on-line b-tagging in fully hadronic 
events were developed specifically for this kind of analysis. Limits on the cross 
section times branching fractions are derived model independently. The result 
was interpreted in the Minimal Supersymmetric Model (MSSM). Upper limits at 
ton/3 are derived as a function of M^, the mass of the Higgs boson A. In the 
analysis presented here data taken in 2011 at 7 TeV are used corresponding to a 
total luminosity of L = 2.7 — 4.0/fe~^ for low and medium Higgs boson masses, 
respectively. 

1 Introduction 

A search is performed for a Higgs boson produced in association with additional 
b-quarks. In the Standard Model (SM) the cross section is very low and due to the 
large background dominated by QCD processes with three b-quark or two b-quark 
plus c,udsg-quark final states the signal in the current data is not in reach. 

In the Minimal Super-symmetric Model (MSSM) [l|[2], there are 3 neutral 
Higgs bosons, one CP-odd (A) and two CP-even {h, H). Their couplings to d-type 
fermions are enhanced by a factor tan(3. The Higgs bosons A and h{H) degenerate 
in the mass. Hence, the cross section of Higgs boson production in association with 
b quarks increases by a factor 2 x tail? (3. 

The MSSM Higgs bosons production is dominated at hadron colliders by the 
sub-process gg — > &(&)$ + X with only a small contribution from qq — > hh^ + X. 
The branching fraction of the Higgs boson decay to b-quarks are in the range 
[70%, 92%] [3], and it remains large even at large masses of the Higgs boson A. 

The dominant background arises from heavy flavor multi-jet QCD processes. 
Other background contributions like Zhh, Wbb, tt have low rates. 

Data collected with the CMS detector ^ in proton-proton collisions at a centre- 
of-mass energy of 7 TeV at the LHC, corresponding to an integrated luminosity of 
2.7 — 4.0/6"^, were used in this search. Specialized triggers have been developed, 
using algorithms to tag b-quark jets from displaced vertices of B-meson decays. 

2 Trigger and event selection 

At LHC multi-jet events are produced copiously via higher order QCD processes. 
To trigger events with jets originating from the decay of heavy particles thresholds 
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on transverse momenta (pr) of the jets are imposed in the first level trigger. In the 
Higgs boson low-mass search, performed in the data taken in the first period, the 
thresholds were idGeV/c for the leading jet and 38 GeV/c for the second leading 
jet. In the second half of 2011, due to the higher luminosity, the threshold for 
the leading jet was 60 GeV/c, and for the second leading jet 53 GeV/c. This data 
were used in the search for the medium-mass Higgs boson. On-line b-tagging was 
performed with TCHE algorithm [H] using the track with 2"'' highest significance, 
and at least two jets must be b-tagged. Tracks were clustered into on-line jets [6] 
with a cone radius i? = 0.5. 

The jets were reconstructed using the particle-flow algorithm I^{9] and a stan- 



dard jet energy calibration was applied 10 . Then, jet energies were corrected for 



the contributions of pile-up (PU) interactions 11 . The combined secondary vertex 



(CSV) algorithm 12 was applied for the offline identification of b-jets. Secondary 



vertices (SV) were reconstructed using the Adaptive Vertex Finder 13 . It combines 



several topological and kinematic variables of the SV and track impact parameters 
into a ratio of likelihoods. The ratio was required to be larger or equal to 0.89 (tight 
operation point, CSVT) to tag the jet. Mistagged light flavor jets are suppressed 
to a level of ^ 0.2% with a b-jet identification efficiency of about 55% [l4] . 

A variable (EventBta) was created exploiting the SV mass spectra obtained for 
jets of different flavors. EventBtag combines the SV masses of the three leading 
jets in one quantity which is small for c- and light-jets and large for b-jets. 

The three leading jets have to satisfy the following criteria in the search of the 
Higgs bosons with low (medium) masses: prijeti) > 46(60), PT{jet2) > 38(53), 
PrUeta) > 20 GeV/c; \ri{jets)\ < 2.2; at least three b-tagged jets with CSVT are 
required. 



3 Background Modeling 

The signature of the process pp hH + V ^ 35 + V is at least three b-jets 
which are identifled with the b-tagging algorithm. The dominant background is 
due to multijet QCD processes mainly composed of events with three b-jets and 
events with two b-jets plus a charm- or udsg(light)-mistagged jet. Contributions 
from other processes like tt, Z-|-jets were estimated to be negligible. 

We have used a data sample with two b-tagged jets to estimate the back- 
ground in data with three b-tagged jets. The double-btag data events were properly 
weighted assuming the flavor (b,c,udsg) of the not-tagged jet to obtain the shape 
of QCD background. Then a 2D template spanned by the di-jet mass of first 
two leading jets, M12 and EventBtag was created by filling with the events of the 
given category. Finally, five templates were used. Their projections in M12 and 
EventBtag are shown in Figure [ij Two additional corrections are applied to re- 
move contamination from non-bb events in the double-btag sample and to account 



for differences between offline and on-line b-tag efficiencies 15 



Signal templates were obtained for each mass point from samples of Higgs bosons 



trigger efficiencies 15 were applied. 



events simulated using PYTHIA 16 . Pileup reweighting and Level- l/Level-2 
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Figure 1. M12 (a) and EventBTag (b) projections for the five background templates used in the 
search for a low-mass Higgs boson. 
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Figure 2. (a) Results of the 'background-only' fit in the triple-btag sample in the projection on 
the di-jet mass, (b) Results of the fit in the di-jet mass projection including a signal template for 
a Higgs boson with a mass of 200 GeV/c^ in the triple-btag data for the medium-mass scenario. 



4 Results 

The signal and background yields were obtained from a binned unconstrained least- 
squares fit of '2D' templates to the triple-btag data distributions in the M12 and 



EventBtag space 15 . The statistical uncertainties of the templates are taken into 
account. 

The fit to the data distribution was performed using only background templates. 
The result for the low-mass scenario of the Higgs bosons search is shown in Figure 

The data distribution is well described by the fit, with a x^/^dof ~ 1- The 
data points are within the uncertainty of the fit, indicated by the cyan area. The 
same 'background-only-fit' for the medium-mass scenario also perfectly describes 
the data distribution. In both cases events with at least three b jets constitute the 
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Figure 3. (a) The observed 95% C.L. upper limits on the cross section multiplied by the branching 
fraction as a function of M^. (b) Observed and median expected 95% C.L. upper limits on tanfi 
as a function of Ma in the mJJ""^ scenario for /i = +200 GeV/c^ 



largest contribution to to the background. For the low-mass Higgs bosons search 
it amounts to about 74%, in good agreement with the MC background estimation 
of 72%. 

The fits were repeated including signal templates for several Higgs boson masses 
between 90 to 350 GeV jc? . As an example, the result of the fit assuming a Higgs 
boson with a mass of 200 GeV/c^ is shown in Figure [2jb). The values of the 
signal fractions obtained in the fits using all Higgs boson mass hypotheses are used 
to derive signal cross sections multiplied with the branching fraction of the Higgs 
boson in a pair of b-quarks as shown in Figure |3][a) . 

Upper limits on the cross section multiplied by the branching fraction, a{pp — ^ 
6$) X ^ hi), were calculated in the [90, 350] GeF/c^ range for Ma- The 

modified frequentist CLs method 17 is applied using the RooStats package [Ts) . 
As test statistics for the CLs method the profile likelihood ratio distribution was 
used. Systematic uncertainties were taken into account as nuisance parameters of 
the likelihood fit to the data. Normalizations of the templates were unconstrained 
in the fit. The observed and the median expected 95% C.L. limits as functions of 
the Higgs boson mass are shown in Figure [3]jb). The observed limit is well within 
the expected 2a band. 

The model independent limits on a{pp — >■ 6$) x BR{^ — )■ bb) were translated 



into tan[3 limits within the m™"^^ scenario 19 , where the parameter fi is set to 
+200 and -200CeV/c^. 

The expected cross section and branching fraction, in the MSSM framework, 



are calculated by BBH@NNLO j20j and FEY N HIGGS [2l]- [24], respectively. 

The observed and expected median 95% C.L. upper limits on tan(3 as a function 
of Ma using /i = +200GeF/c^ are shown in Figure pi b). The upper limits on tanj3 
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for = —200 GeV/c^ case were compared with the recent results from DO [25^ and 
CDF [26^ experiments at the Tevatron coUider. The limits obtained in the current 
analysis are significantly lower than the Tevatron results. 
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